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A  COMPUTER  PROGRAM  FOR  ESTIMATION  OF  PARAMETERS  FOR  THE  WEIBULL 
INTENSITY  FUNCTION  AND  FOR  THE  CRAMER-VON  MISES 
GOODNESS  OF  FIT  TEST 

1.  INTRODUCTION 


The  Weibull  intensity  function 

U(X)  =  XBX6"1  (1 .1 ) 

X>0,  3>0,  X>0,  is  frequently  used  as  a  model  for  the  determination  of 
reliability  growth  and  wear-out  characteristics  for  a  wide  variety  of 
complex,  repairable  systems.  The  failure  rates  of  military  equipment 
such  as  vehicles,  aircraft,  guided  missiles,  electronic  computer  systems, 
and  ammunition  are  being  evaluated  using  this  model. 


Formulas  have  been  developed  by  the  Army  Materiel  Systems  Analysis 
Activity  (AMSAA)  for  the  maximum  likelihood  estimation  of  the  unknown 
parameters  x  and  6,  based  upon  sample  data.  These  estimation  formulas 
for  the  Weibull  process,  found  in  Crow  (1975),  can  be  stated  as 


Km 

INi 

i  =1 


*  Z  [T?21n(Ti2)-TJi1nCT1i)] 


K  N. 

I  V 

1=1  j*l 


(1.2) 
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where: 


(1.3) 


K  is  the  number  of  systems  under  study; 

N.  is  the  total  number  of  failures  (or  occurrences  of  an  event 
under  study,  such  as  unscheduled  maintenance  actions;  etc.) 

for  the  iin  system; 

T.,  is  the  starting  time  of  the  period  of  continuous  observation 

1  I  it. 

of  the  iin  system; 

T.?  is  the  ending  time  of  the  period  of  continuous  observation 

of  the  iLn  system; 
t  h 

X..  is  the  jLn  time  of  occurrence  of  the  failure  (or  event),  for 

1 J  th 

the  i  system; 

In  is  the  natural  logarithm,  and  0-ln(0)  is  defined  to  be  0. 


i  I  •* 


f  f  if 
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To  expedite  the  computation  of  these  estimates,  AMSAA  developed  a 

FORTRAN  computer  program  to  calculate  3  and  X.  That  program,  which  was 

documented  in  Belbot  (1974),  was  successfully  employed  by  the  U  S  Army 

Materiel  Development  and  Readiness  Command  (DARCOM),  various  subordinate 

commands  and  several  project  managers'  offices,  as  well  as  by  AMSAA. 

After  the  parameters  were  determined  by  that  program,  a  goodness  of  fit 

test  was  frequently  used  to  test  statistically  the  hypothesis  that  the 

failure  times  of  the  systems  being  analyzed  followed  a  nonhomogeneous 

Poisson  process  with  Weibull  intensity  function  (see  Crow[  1975]).  The 
* 

modified  Cramer-Von  Mises  goodness  of  fit  statistic  was  computed,  either 
by  hand  or  by  a  separate  computer  program. 

Obviously,  the  consolidation  of  an  automated  goodness  of  fit  test  with 

the  computer  routine  which  estimates  the  parameters  3  and  X,  would  increase 
efficiency  and  accuracy.  Either  the  manual  calculations  or  the  use  of  a 
separate  program  would  be  eliminated.  To  accomplish  properly  this  consol ida 
tion,  the  estimation  procedure  was  subordinated  to  a  new  main  program  which 
also  controls  the  input  of  data  and  the  execution  of  the  goodness  of  fit 
test.  Because  of  the  radical  nature  of  this  redesign,  it  was  appropriate  to 
incorporate  other  new  features  at  the  same  time.  Principal  among  these  new 
features  are  a  simplified  input  procedure  and  dynamic  data  storage  alloca¬ 
tion.  The  resulting  computer  program  is  easier  to  use  and  provides  more 
information  that  its  antecedent  program.  This  report  will  explain  the 
structure  and  the  use  of  this  new  program. 

2.  COMPUTING  PROCEDURE 

2.1  Estimation  of  Parameters 


Since  the  formulas  (1.2)  and  (1.3)  do  not,  in  general,  yield 

A  A 

&  and  X  in  closed  form,  an  iterative  technique  is  required.  Formula 
(1.2)  may  be  recast  as 


K  N. 

I  l  ln(X  ) 
i=l  j=l  1J 
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(2.1) 


by  replacing  X  by  its  equivalent  expression  from  equation 
execution  of  a  few  simple  algebraic  operations.  Equation 

consists  of  a  constant  with  regard  to  3,  minus  a  function 
0,  or  simply 


(1.3),  and  by 
(2.1)  now 

/V 

of  3,  yielding 


A  -  D(3)  =  0 


(2.2) 


2 


The  correct  value  of  3  will  satisfy  equation  (2.2)  and  can  be  used  to 

A 

calculate  the  corresponding  value  of 

/V  I 

The  solution  for  3  is( iteratively  determined  in  the  following  way. 

For  an  initial  estimate  3  which  is  assumed  to  be  larger  than  the  true  3, 

the  expression  A-D(3  )  is  evaluated.  For  all  values  of  3  larger  than 

the  true  3,  the  subtraction  yields  a  negative  result.  After  each  negative 

result,  3  is  reduced  by  the  initial  step  size  of  1,  and  A-D(3  )  is  again 
evaluated. 

(  When  a  positive  number  results  from  the  subtraction,  indicating  that 

3  is  too  small,  the  step  size  is  decreased  to  0.10  of  the  present  step 

size,(the  previous  value  of  3  which  gave  a  negative  result  for 

A-D ( 3  ),  is  adjusted  by  the  new  step  size  and  the  evaluation  process 
begins  again. 

The  iteration  procedure  continues,  adjusting  3  by  the  new  step  sizes, 
until  the  left  side  of  equation  (2.2)  is  within  a  specified  tolerance  E  of 

0.  *  is  then  calculated  based  on  3,  using  factors  already  computed  in 

finding  3.  This  procedure  is  summarized  by  the  state  diagram  (Figure  2.1). 

2.2  Goodness  of  Fit  Test 


The  Cramer-Von  Mises  Goodness  of  Fit  Test  is  appropriate  whenever 
the  starting  time  of  each  system  is  equal  to  0.  To  perform  this  test, 
the  program  first  transforms  the  failure  times.  For  time  truncated  test¬ 
ing,  the  failure  times  for  each  system  are  divided  by  the  ending  time  of 
the  test  period  for  that  system.  In  failure  truncated  testing,  for  every 
system,  all  the  failure  times  except  the  last,  are  divided  by  the  last 

failure  time.  The  last  failure  time  is  thereafter  excluded  from  the 

calculations  and  the  number  of  transformed  failures  is  one  less  than  the 
original  number  of  failures  for  each  such  system.  All  the  transformed 
failure  times  are  then  sorted  into  increasing  order. 

Next,  the  unbiased  estimate  6  of  the  estimated  shape  parameter 
&,  is  calculated  using  the  relation: 

6  =  - ^ - - -  (2.3) 

K  k  L . 

I  L..ln(T.2)  -  l  rin(X..) 

i=l  1  u  i=l  j=l  1J 

( 

where: 

j_  i_ 

!N - — 1  if  the  data  for  the  i  system  are  failure  truncated 
-  .th 

N.j  if  the  data  for  the  i  system  are  time  truncated 
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Figure  2.1  State  Diagram 


and  M  =  total  number  of  transformed  failure  times.  (Note  that 

M  =  i2ili •> 

2 

Finally,  the  Cramer-Von  Mises  statistic  CM  is  computed  by  the 

formula: 

i  ■  m  ♦  j,  <z?  -  <2-<> 

where  the  Z.  are  the  transformed  failure  times.  An  explanation  of  this 
i  2 

test  and  a  table  of  critical  values  of  CM  may  be  found  in  Crow  (1975). 


Next  page  is  blank. 
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3.  DESCRIPTION  OF  PROGRAM 


3.1  Major  Features 

The  program,  which  is  listed  in  Appendix  A,  has  some  important 
features.  First,  the  program  is  written  in  American  National  Standards 
Institute  (ANSI)  FORTRAN  X3. 9-1966,  and  should  therefore  execute  on  any 
computer  having  a  compiler  for  this  language.  Secondly,  the  amount  of 
storage  required  to  use  this  program  should  not  cause  difficulties  since 
all  data  arrays  are  dynamically  allocated  under  control  of  the  main 
program  (see  Chung-Phil lips,  et  al.,  1975  ). 

Finally,  and  perhaps  most  importantly  for  the  user,  this  program 
uses  free-field  input,  that  is,  no  specific  format  is  required  for  the 
input  data.  The  time  and  the  effort  regularly  expended  in  preparing 
data  for  input,  are  greatly  reduced  because  of  this  feature.  Moreover, 
the  misalignment  of  data  fields  to  formats,  a  frequent  source  of  errors 
in  using  many  computer  programs,  is  eliminated  entirely.  While  no  input 
scheme  can  be  regarded  as  foolproof,  free-field  input  is  much  more 
flexible  than  fixed-field  format  specifications. 

3.2  Overall  Characteristics 

All  calculations  in  the  program  are  made  in  double  precision 
mode.  Experience  has  shown  that  the  use  of  single  precision  variables 
for  these  calculations  often  results  in  significant  discrepancies  in  the 
estimates  of  the  parameters  due  to  errors  accumulated  during  the  iterative 
process. 


The  modular  structure  of  the  program  (see  Figure  3.1)  reflects 
organization  by  functional  purpose.  The  input  of  data,  certain  inter¬ 
mediate  calculations,  the  estimation  of  parameters,  and  the  goodness  of 
fit  test  are  each  performed  by  an  independent  module.  Major  subprograms 
print  their  results  as  the  values  become  available.  Subroutines  which 
detect  errors,  print  diagnostic  messages  naming  the  detecting  routine 
and  briefly  stating  the  difficulty,  and  then  attempt  to  continue  process¬ 
ing  when  possible.  Independence  of  the  subroutines  is  maintained  by 
restricting  communication  between  individual  subprograms  to  the  passage 
of  formal  parameters  in  argument  lists.  No  COMMON  statements  are  used. 

3.3  Specific  Details  of  Routines 

In  addition  to  controlling  the  major  modules,  the  main  program 
also  allocates  storage  for  data  arrays,  as  stated  in  Section  3.1.  The 
allocation  is  based  upon  the  maximum  number  of  failures,  NFAIL,  and  the 
maximum  number  of  systems,  NSYS.  The  master  data  array,  BLOCK,  has  a 
length  of  NTOTAL,  equal  to  1.5  times  the  value  of  NFAIL  plus  six  times 
the  value  of  NSYS.  To  redimension  the  entire  program,  one  merely  adjusts 
the  values  of  NFAIL,  NSYS  and  NTOTAL  in  the  DATA  statements  at  the  beginning 
of  the  main  routine  and  changes  the  size  of  the  array  BLOCK,  also  found 
at  the  beginning  of  the  main  program,  to  equal  the  new  value  of  NTOTAL. 
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Figure  3.1  Module  Organization  Chart. 


By  these  actions,  all  data  arrays  in  all  subroutines  will  be  properly 
resized.  Since  the  program  size  excluding  data  arrays,  is  less  than 
5,000  words,  computer  memory  requirements  can  be  scaled  to  problem  size 
through  use  of  this  feature. 

Also  found  in  the  DATA  statements  at  the  beginning  of  the  main 
routine  are  the  unit  number  for  input,  IUNIT,  and  the  logical  switch  ECHO 
which  controls  the  printing  of  the  input  data.  Just  as  the  storage 
allocation  values,  these  values  may  be  changed  as  needed. 

''  > . 

The  first  major  module  is  for  the  input  of  data.  The  INDATA 
module  reads  the  beginning  and  the  ending  times  and  the  failure  times 
from  the  input  unit,  IUNIT.  If  the  logical  variable  ECHO  is  true,  the 
submodule  REPEAT  will  print  the  input,  using  the  LOCIP  subroutine  to 
isolate  in  storage  the  failures  for  each  system  and  the  PRNTX  subroutines 
to  print  them.  The  failure  times  are  stored  in  a  linear  array  called  X. 
Parallel  to  X,  an  integer  array  named  MARKER  identifies  the  failures 
with  the  respective  systems,  such  that  if  failure  X(I)  occurred  on 
system  J,  then  MARKER(I)  equals  J. 

The  free-f ield  reading  of  data  is  performed  by  the  INREAL  sub- 
module.  This  submodule,  consisting  of  the  routines  INREAL,  CHAR,  INSYML, 
MATCH  and  UNSGND,  is  a  translation  from  ALGOL  into  FORTRAN  of  Algorithm 
239  of  the  Association  for  Computing  Machinery  (see  McKeeman  [1964]). 
Specific  details  concerning  the  input  arrangement  are  given  in  Section  4. 

The  INDATA  subroutine  also  sets  three  logical  variables  depend¬ 
ing  on  the  input.  If  data  errors  are  encountered,  the  variable  NOGOOD  is 
made  true.  If  all  systems  start  at  0.0,  then  the  goodness  of  fit  test 
will  be  appropriate  and  so  the  variable  GOF  is  set  to  true.  Lastly,  when 
the  end  of  input  is  reached,  the  logical  variable  HALT  is  returned  as  true 

The  second  principal  module,  LOGS,  calculates  the  logarithms  of 
the  beginning  and  the  ending  times,  and  the  sum  of  the  logarithms  of  the 
failure  times.  For  computational  purposes,  a  beginning  time  of  0.0  is 
defined  to  have  a  logarithm  of  0.0  instead  of  infinity.  Note  that 
failure  times  of  0.0  are  not  valid  for  this  model  and  are  flagged  as 
errors  by  the  preceding  INDATA  module. 

The  next  major  module,  MLEW,  computes  the  maximum  likelihood 
estimates  of  the  parameters  of  the  Weibull  intensity  function  for  the 
given  data,  using  the  formulas  discussed  in  Section  2.1.  If  unsuccessful, 
this  module  will  report  one  of  three  possible  error  conditions.  The 
first  error  message  "BETA  LESS  THAN  0.00000001"  indicates  that  the 

data  should  be  rechecked  for  the  reasonableness  of  a  very  small  p.  The 
second  error  message  is  "INITIAL  ESTIMATE  OF  BETA  IS  TOO  SMALL."  Since 
the  initial  estimate  of  BETA  is  set  to  10  at  the  start  of  the  MLEW  sub¬ 
program,  this  message  indicates  some  pecularity  of  the  data.  (In 

general,  0<B<10.)  The  third  message,  "STEP-SIZE  HAS  BECOME  INSIGNIFICANT 
BETA  NOT  RESOLVABLE,"  indicates  that  the  module  has  gone  as  far  as 
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possible  trying  to  meet  the  tolerance  set  for  the  difference  A-D(3  ). 

This  tolerance,  EPSILN,  may  be  enlarged  by  changing  the  assignment  state¬ 
ment  also  located  at  the  beginning  of  the  MLEW  subprogram. 

The  last  major  module,  CVM,  performs  the  Cramer-Von  Mises  goodness 
of  fit  test,  as  described  in  Section  2.2.  The  failure  times  for  each 
system  are  located  in  storage  using  the  LOCIP  subroutine,  and  examined  to 
determine  if  the  testing  was  time  truncated  or  failure  truncated.  The 
failure  times  are  then  transformed  by  the  TRANSF  subroutine  and  sorted  by 
the  SORTX  subroutine.  (SORTX  is  a  modification  of  an  utility  subprogram 
described  in  Campbell,  et  al . ,  [1970].)  The  unbiased  estimate  of  3,  UNBETA, 
is  calculated  next,  as  explained  earlier.  The  last  phase  depends  on  the 
system  starting  times.  If  any  starting  time  is  non-zero,  the  module 

terminates  with  a  message  stating  that  the  Cramer-Von  Mises  goodness  of 
fit  test  is  not  appropriate.  Otherwise,  the  goodness  of  fit  statistic, 

CM2,  is  computed  and  printed. 
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4.  INPUT  REQUIREMENTS 


As  stated  previously,  the  input  for  this  program  is  free-field.  The 
only  requirement  regarding  spacing  is  that  at  least  one  blank  column 
separate  adjacent  values.  The  values  must  not  run  together.  This  means 
that  the  program  generally  takes  the  same  view  of  the  data  that  a  person 
would,  namely,  that  each  cluster  of  numeric  characters  constitutes  one 
data  value.  The  only  exception  to  this  rule  occurs  at  the  boundaries  of 
records.  Since  the  input  is  treated  as  a  continuous  stream,  a  string  of 
characters  beginning  in  the  first  column  of  a  record,  is  considered  a 
continuation  of  the  string  of  characters  ending  in  the  last  column  of  the 
previous  record,  if  any.  Record  boundaries  are  not  delimiters;  blanks 
are  the  only  delimiters. 

The  data  required  for  this  program  consist  of  the  beginning  and  ending 
times  of  the  test  period  for  each  system,  and  the  failure  times  for  each 

system.  The  arrangement  of  the  input,  (which  is  also  stated  in  the 

comments  of  the  INDATA  subroutine),  is  by  system.  The  first  data  value 
is  the  beginning  time  of  the  first  system.  The  second  data  value  is 
the  ending  time  of  the  first  system.  Next  is  the  failure  times  for  the 
first  system,  followed  by  a  negative  value  to  mark  the  end  of  the  first 
system.  The  same  pattern,  beginning  time,  ending  time,  failure  times, 

and  negative  trailer,  is  repeated  for  each  subsequent  system  in  the  first 

data  case. 

Another  negative  value  (making  two  in  a  row),  signals  the  end  of 
input  of  the  current  data  case,  and  the  beginning  of  the  computational 
procedures.  The  same  arrangement  may  be  repeated  for  as  many  cases  as 
desired  per  program  run.  When  the  input  routine  encounters  a  negative 
value  after  completing  a  case,  (that  is,  the  third  negative  value  in  a 
row),  the  end  of  the  program  run  is  indicated. 

Thus,  as  a  simple  example,  if  one  desired  to  use  this  program  for 
one  run  consisting  of  one  case  wherein  one  system  experienced  seven 
failures,  the  input  data  would  be:  the  beginning  time,  the  ending 
time,  the  seven  failure  times,  and  three  negative  values.  To  demonstrate 
the  latitude  of  the  input  requirements,  the  data  for  a  number  of  test 
cases  are  shown  in  Appendix  B.  Notice  that  any  negative  value  is  accept¬ 
able  as  a  trailer  and  that  data  values  may  be  entered  with  or  without 
decimal  points.  Although  not  shown  in  the  examples,  data  values  may 
also  be  in  exponential  form,  that  is,  containing  ' E , 1  '+'  or 


Next  page  is  blank. 


11 


5.  TRANSFERABILITY  AND  MODIFICATION 


Since  this  program  is  written  in  standard  FORTRAN,  transfer  to  other 
computer  systems  should  be  straightforward.  To  assist  in  the  transfer 
process,  Appendix  C  contains  the  output  produced  by  the  program  for  the 
input  shown  in  Appendix  B.  This  output  was  generated  on  a  Control  Data 
Corporation  (CDC)  Cyber  76  Computer,  using  the  program  exactly  as  listed 
in  Appendix  A.  (Note  that  non-standard  PROGRAM  statement  required  by 
the  CDC  Cyber.) 

The  input  for  these  test  cases  came  from  records  of  eighty  characters 
each.  If  the  input  record  length  is  other  than  eighty,  two  changes  may 
be  required.  The  value  of  the  variable  LENGTH  and,  if  necessary,  the 
dimensioned  size  of  the  array  BUFFER,  should  be  adjusted  in  the  INSYML 
subroutine. 

Alternatively,  one  could  replace  the  entire  INREAL  submodule.  Although 
these  routines  were  written  to  be  fully  transportable,  running  time  might 
be  saved  by  using  the  system  defined  free-field  reading  capability  of  any 
computer  having  such  a  feature.  As  an  example,  a  substitute  for  the 
INREAL  submodule,  suitable  for  the  CDC  Cyber  76,  is  shown  in  Figure  5.1. 
Such  substitutes,  however,  are  system  dependent  and  not  readily  trans¬ 
ferable. 
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FIGURE  5.1  SUBSTITUTE  FOR  INREAL  SUBMODULE 


SUBROUTINE  INREAL  (IU,  X) 

DOUBLE  PRECISION  X 

READ  (IU,  *)  X 

RETURN 

END 
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BEGINNING  AND  ENDING  TIMES  FOR  FIRST  SYSTEM,  FAILURE 

TIMES  FOR  FIRST  SYSTEM,  NEGATIVE  VALUE  AS  TRAILER. 
BEGINNING  AND  ENDING  TIMES  FOR  SECOND  SYSTEM,  FAILURE 

TIMES  FOR  SECOND  SYSTEM,  NEGATIVE  VALUE  FOR  TRAILER. 
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INPUT 

DATA 

e.o 

-999. 

2C0.  0 

6.3 

6.6 

10.2  23.5 

23.8 

26.6 

76.0 

77.1  92.1 

197 

.2 

w  •  \  ’ 

-999 

2  •»  C  •  > 

5.1 
12  . 

5.6 

1  161 

18.6  19.5 

.8  18,i.6 

26.2 

19i0, 8 

26.7 

65.1 

65. e  75.7 

79. 

7 

98.6 

*  • 

-9Q9 

2*".  .  > 

8.6 

195. 

32.5 

6 

66.7  68.6 

50.6 

73.6 

98.  7 

112.2  129. 

8 

136. 0 

-999 

. 

-999. 

197.2 

6.3 

6.6 

10. 2  23*5 

23,8 

26,6 

76.0 

197.2  77.1 

92 

•  1 

C.r.  • 

—  999 

190.  8 

€•  1 
120# 

5.6 

1  161 

18.6  19.5 

.8  180.6 

26.2 

19C.6 

26,7 

65.1 

65.8  75,7 

79. 

7 

98.6 

19?.  8 

195.8 

-999 

-999 

0.4 

32.5 

66.7  68.6 

50.6 

73.6 

98.7 

112.2  129. 

6 

136 

. 

r  i(.  -3<? 

1? 

-99  U 

2.  6 

12  17  -99  0 

20  -99 

0  1C 

-99 

-99 


liK.J  32.  66.  56.  75.  95.  -99.  -99 

12.2  6:‘,3  10. C  15.5  16.8  857.9  -99  -99 

3256.3  .7  3.7  13.2  17.6  56.5  99.2 
176.*  191.6  262.8 

355.2  686.  3  692?. 5  513.3  558.6  678.1  688. C 

1.  ?«.l  H36.6  1136.1  1178.9  1259.7  1297.9 

1629. 8  1702.3  1928.9  2072.3  2525.2  2928.5 

3256.3 


112.2  120.9  151.0  163. C 

786.9  8 87.0  1010.7 

1619.7  1571.7 

3016.6  3181.0 


-99.  -99. 

C>.2'  79, a  >,c  20  65.  62  86  79  68  -8  -7 
656.  789.  656  656.1  657  658  689  -99  -99 

0.0  800. 0  65.  65'0»  667.  677.  686.  692. 
v.12  566.0  65.0  78.  89.  99.1  -99  -99 

.  lv .  2.  3.  6,  5.  6.  8.  -99.  )%  30.  -99.  C.  50.  61.  62.  63.  66.  -99 

0.  *w.  51.  52.  53.  -99.  -99. 

-99 
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WEIBULL  INTENSITY  MODEL 


PARAMETER  ESTIMATION  AND  GOODNESS  OF  FIT  TEST 


VERSION  OF  12  NOV.  1979 


37 


DATA  INPUT  PHASE 


TOTAL  NUMBER  OF  FAILURES  -  36 

TOTAL  NUMBER  OF  SYSTEMS  ■  3 


SYSTEM  STARTING  AND  ENDING  TIMES. 


FAILURES. 


2  *0000+02 

4 • 300D+00 

4.4000*00 

1.0200*01 

2.3500*01 

2.3800*01 

2  *6400  +  01 

7. 4000*01 

7.7100+01 

9.2100*01 

1.9720*02 

2  *0000+02 

1.0000-01 

5.6000*00 

1 . 860D  +  01 

1 *9500+01 

2  *  420D+01 

2.6700+01 

4.5100*01 

4.580D+01 

7.570D+01 

7.970D+01 

9.8600+01 

1.2010*02 

1.618D+02 

1.606D+02 

1.908D+02 

2*0000+02 

0. 4000*00 

3.2500*01 

4.4700*01 

4.8400*01 

5 .0600*01 

7.3600*01 

1.958D+02 

9.8700*01 

1.1220*02 

1.2980*02 

1.3600*02 

33 


ESTIMATION  OF  THE  PARAMETERS  OF  THE  WEIBULL  INTENSITY  FUNCTION 


CONSTANT  NOT  INVOLVING  8ETAI  A  •  3.67318960+00 


ESTIMATED  BETA 

1.000000000D+01 
9. 00000000 OD +00 
8.0000000000+00 
7 . 00000000 OD +00 
6.0000000000+00 
5 • 000000000 D +00 
4 .0000000 OOD+OO 
3 • 000 000000 D +00 
2. 000000000 D+00 
l.OOOOOOOOOD+OO 
9.000000000D-01 
8.000000000D-01 
7.0000000000-01 
6. OOOOOOOOOD-Ol 
6 .900000 000 D-01 
6.800000000D-01 
6 .700000000 D-01 
6.600000000D-01 
6.500000000D-01 
6 .400000000 D-01 
6 .3 00000 000 D-01 
6 .200 000000 D-01 
6 . 1 00000 OOOD— 01 
6.1900000000-01 
6. 180000000 D-01 
6.170000000D-01 
6.160000000D-01 
6.150000000D-01 
6.159000 000D-01 
6. 158000000 D-01 
6  * 157000000D-01 
6.1560000000-01 
6. 155000000 D— 01 
6.1540000000-01 
6.153000000D-01 
6.1539 00 OOOD -01 
6.153800000D— 01 
6.  153700000D-01 
6. 153 600000 D-01 
6*153500 OOOD -01 
6.153400000D-01 


FUNCTION  D(8») 

5 • 1 983 1736 7D+00 
5 . 187206255  D+00 
5.17331 7367D +00 
5 . 15 5 460 22 4 D+00 
5.1316507000+00 
5.0983173670+00 
5. 048317 36 7D +00 
4 .96498 40 33 D+00 
4.798317 36 7D+00 
4 .298317 36 7D +00 
4.187206255D+00 
4 .048 31 736 7D +00 
3 .8 69745 93 8 D+00 
3.631650700D+00 
3.849042004D+00 
3.827729131D+00 
3.805780053D+00 
3.783165851 D+00 
3 .75985  5 828 D+00 
3.735817367D+00 
3 ,7110 15779 D+00 
3 .685414141 D+00 
3 .658973 10 4 D+00 
3.682808481D+00 
3.680194389D+00 
3 .677571 82 4D +00 
3. 674940 743D+00 
3.672301106D+0Q 
3 .6 7467716 5D +00 
3.674413502D+00 
3 .6741497520+00 
3 .673885918 D+00 
3.673621997D+00 
3 .67 33 5 7 99 ID +00 
3 .673093  89 8D+00 
3 .6 73 3 31 58 5 D+00 
3.673305179D+00 
3 .673278772 D+00 
3.673252364D+00 
3.673225955D+00 
3 .673 19 9 54 6 D+00 


A  -  D(8*  ) 

-1. 5251277180+00 
-1. 5140166070+00 
-1.5001277180+00 
-1.4822 70 575D+00 
-1.4584610510+00 
-1. 4251277180+00 
-1.375127718D+00 
-1 .2917943840+00 
-1.125127718D+00 
-6. 2512771780-01 
-5 . 140 166067 D-01 
-3.7512771780-01 
-1.965562892 D-01 
4. 153894890 D-02 
-1. 7585235540-01 
-1.545394825D-01 
-1 .3259040430-01 
-1.0997620260-01 
-8 .666617931 D-02 
-6 .262 7 71 777 D-02 
-3 .7826 13 047 D-02 
-1.2 22449196 D-02 
1.42165 445 3 D-02 
-9.6188 32469 D-03 
-7. 0047404 22 D-03 
-4. 3821748 19 D-03 
-1. 751094 39 2D-03 
8. 885 423943 D-0 4 
-1 .4 87 5 164 3 7 D-03 
-1 .22385287 7D-0 3 
-9 .6010367050-04 
-6 . 962 687754 D-0 4 
-4.323481501D-04 
-1.683417527D-04 
9.575045862D-05 
-1.41 936393 8 D-0 4 
— 1. 155301769D-04 
-8 .912310164 D-0 5 
-6 .271 51681 6 D-0 5 
-3.630637637D-05 
-9.8967262360-06 


ESTIMATED  LAMBDA 

1 .1718750000-22 
2 .3437500000-20 
4. 68 7 500000 D- 18 
9.3 75 000000 D- 16 
1 .8750000000-13 
3.750000000D-11 

7 .5000000000- 09 

1.5000000000- 06 
3.0000000000-04 
6.000000000D-02 
1.0191 87879 D-01 
1.731239887 D-01 
2 .9407645140-01 
4.995319244D-01 
3 .100 777107 D-01 
3.269496291D-01 
3 .4473958060-01 
3.634975173D-01 
3 .8 3276 10 89D -01 
4.0413069 14 D-01 
4.261204223D-01 
4.4930 64 454 D-01 
4.7375406420-01 
4 .516933312 D-01 
4 . 540 92 89 71 D-01 
4.565052103D-01 
4. 589303387 D-01 
4.613683502D-01 
4.5917355900-01 
4. 5941 690 82 D-01 
4 .596603 863 D-01 
4. 5 99 03 99 35 D-01 
4.601477298 D-01 
4. 60 391 59 5 3 D-01 
4 .60 63 55900 D-01 
4. 6041 598 89 D-01 
4. 604403 839 D-01 
4. 60464780 ID -01 
4. 60 4 8 91 7 76 D-01 
4.605135765D-01 
4.605379766D-01 
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THE  FINAL  ESTIMATE  OF  BETA  «  6.1534000D-01. 

THE  FINAL  ESTIMATE  OF  LAMBDA  ■  4 .6053798D-01 

CONVERGENCE  TO  9.8967262D-06 

WHICH  IS  LESS  THAN  EPSILON  -  1 .00000000-05 

THE  FINAL  STEP  SIZE  IS  1 .0000000D-05 
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CRAMER  -  VON  MISES  GOODNESS  OF  FIT  TEST 


ESTIMATED  BETA  •  6.1534000D-01 

NUMBER  OF  FAILURES  -  36 

UNBIASED  ESTIMATE  OF  BETA  ■  5 .982435B0-O1 


CRAMER  -  VON  MISES  STATISTIC  ■  6.95308910-02 

REJECT  THE  WEIBULL  INTENSITY  MODEL  IF  THE  STATISTIC 
EXCEEDS  THE  APPROPRIATE  CRITICAL  VALUE  FOR  M  -  36 
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DATA  INPUT  PHASE 


TOTAL  NUMBER  OF  FAILURES  ■  36 

TOTAL  NUMBER  OF  SYSTEMS  -  3 


SYSTEM  STARTING  AND  ENDING  TIMES. 


FAILURES. 


1.972D+02 

4  •  300D+00 

4.400D+00 

1.020D+01 

2.350D+01 

2 . 360D+01 

2.640D+01 

7.400D+01 

1.972D+02 

7.710D+01 

9.210D+01 

1.908D+Q2 

l.OOOD-Ol 

5 .6000+00 

1  •  860D+01 

1.950D+01 

2  •  420D+01 

2.670D+01 

4.510D+01 

4.580D+01 

7. 570D+01 

7.970D+01 

9.860D401 

1  .201 D+02 

1 *6 18D+02 

1. 606D+02 

1.908D+02 

1.958D+02 

1.958D+02 

8  •  ^OOD+OO 

3  *2  50D  +  01 

4.470D+01 

4 •  840D+01 

5  *060D+01 

1  •  360D  +  02 

7.360D+01 

9.870D+01 

1 .122D+02 

1 .298D+02 
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ESTIMATION  OF  THE  PARAMETERS  OF  THE  WEIBULL  INTENSITY  FUNCTION. 

CONSTANT  NOT  INVOLVING  BETA!  A  -  3.67318960+00 

ESTIMATED  BETA  FUNCTION  D(B»)  A  -  D ( B  *  )  ESTIMATED  LAMBDA 


1.000000000D+01 
9.000000000D+00 
8. 000000000 D+00 
7.000000000D+00 
6. 0000 00000 D+00 
5 • 000 000 00 0D+ 00 
4.000000000D+00 
3. 00000000 OD +00 
2 .0000000000+00 
l.OOOOOOOOOD+OO 
9.000000000D-01 
8. 000000000 D-01 
7 .000000 000 D -01 
6 .0000000000-01 
6 • 900000000D— 01 
6.8000000000-01 
6.7000000000-01 
6 . 600000000 D-01 
6. 5000000 00 D-01 
6. 4000000 00 D-01 
6 • 300 000 000 D-0 1 
6 • 200 000 000 D-01 
6.2 90000000 D-01 
6.2800 00 000 D-01 
6 • 270000000 D— 01 
6. 2600000000-01 
6 .2 50 000000 D-01 
6.259000000D-01 
6. 258000000 D-01 
6. 2589 00 0 00 D— 01 
6.258800000D-01 
6. 258700 000 D-01 


5 .1727643540+00 
5.161470246D+00 
5.147396282D+00 
5.129352022D+00 
5,1053533810+00 
5 .071828974 D+00 
5.021635991D+00 
4.938107813D+00 
4,7712444930+00 
4  .271046084D  +  00 
4.159915038D+00 
4, 021 00 6 19 7D +00 
3.842414799D+00 
3.6042995750+00 
3, 821708 86 8 D+00 
3, 8 00 3 9 399 7D +00 
3.778442920D+00 
3.7558267200+00 
3 .73 25 14 69 8 D+00 
3.708474238D+00 
3.683670652D+00 
3.658067015D+00 
3.681146920D+00 
3,6786151520+00 
3.676075308D+00 
3.6735273500+00 
3. 670971 24 OD +00 
3.673272107D+00 
3.673016782D+00 
3.673246578D+00 
3  .67322 1 04 8D +00 
3.673195518D+00 


-1.4995747050+00 
-1.488280598D+00 
-1.474206633D+00 
-1.456162373D+00 
-1.432163732D+00 
-1.3986393260+00 
-1.348446342D+00 
-1.264918165D+00 
-1,0980548440+00 
-5. 978564354 D-01 
-4.8672538910-01 
-3.478165480D-01 
-1. 6922515040-01 
6, 889007354 D- 02 
-1. 485192188D-01 
-1.2720434790-01 
-1.052532 71 6D-01 
-8.263707150D-02 
-5,9325049690-02 
-3, 528458948 D- 02 
-1.0481003340-02 
1 . 51 226341 7 D-0 2 
-7.957271513D-03 
-5, 4 25 50 29 5 7 D-0  3 
-2. 885659220 D-0 3 
-3, 37 70 160 37 D- 04 
2. 21 840 8839 D-0 3 
-8, 24 5 8 081 3 5 D-0 5 
1.728670081D-04 
-5  •  69292435  7 D-0 5 
-3.139959008D-05 
-5.869120840D-06 


1.5274022350-22 
2.9774619750-20 
5.8030872290-18 
1, 130 81393 4D -15 
2.203137001D-13 
4.2915017950-11 
8.3578330940-09 
1.627398353D-06 
3 ,1681740720-04 
6.166495375D-02 
1. 0446 15 82 4D -01 
1.7695951690-01 
2.997715275D-01 
5  .078155206  0-01 
3.159962865D-01 
3.3309918360-01 
3.5112774620-01 
3.7013207410-01 
3, 901649784 D-01 
4.1128212870-01 
4.335422075D-01 
4.5700707330-01 
4,3583342960-01 
4.381367604D-01 
4. 4 045 226 40 D-01 
4.427800047D-01 
4.451200472D-01 
4,4301345430-01 
4 .43 2 470 2 70 D-01 
4. 43036806 ID -01 
4, 430601590 D-01 
4.430835132D-01 
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THE  FINAL  ESTIMATE  OF  BETA  ■  6.25870000-01 

THE  FINAL  ESTIMATE  OF  LAMBDA  -  4.43003510-01 

CONVERGENCE  TO  5.8691208D-06 

WHICH  IS  LESS  THAN  EPSILON  •  1.0000000D-05 

THE  FINAL  STEP  SIZE  IS  1 .OOOOOOOD-Q5 
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CRAMER  -  VON  MISES  GOODNESS  OF  FIT  TEST 


ESTIMATED  BETA  ■  6 .2 587000D-01 

NUMBER  OF  FAILURES  -  36 

UNBIASED  ESTIMATE  OF  BETA  >  5 .57260480-01 


CRAMER  -  VON  MISES  STATISTIC  -  1. 3643898D-01 

REJECT  THE  WEIBULL  INTENSITY  MODEL  IF  THE  STATISTIC 
EXCEEDS  THE  APPROPRIATE  CRITICAL  VALUE  FOR  M  «  33 
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DATA  INPUT  PHASE 


TOTAL  NUMBER  OF  FAILURES  -  3 

TOTAL  NUMBER  OF  SYSTEMS  -  5 


SYSTEM  STARTING  AND  ENDING  TIMES. 


0. 

0. 

0. 

0. 

0. 


1.000D+01 
1 . 500D+01 
2.000D*01 
2.000D+01 
1.000D+01 


FAILURES. 

NO  FAILURES 
NO  FAILURES 

4  .0030 ♦00  1. 2000*01  1.700D+01 
NO  FAILURES 
NO  FAILURES 


ESTIMATION  OF  THE  PARAMETERS  OF  THE  WEI8ULL  INTENSITY  FUNCTION 


CONSTANT  NOT  INVOLVING  8ETA«  A  •  2 .23480 48D+00 


ESTIMATEO  8ETA 

1.000000000D+01 
9.0000000000+00 
8.0000000000+00 
7  .000000 OOOD+OO 
6.0000000000+00 
5.000000000D+00 
4.0000000000+00 
3.000000000D+00 
2  .0000 00 OOOD  +  OO 
l.OOOOOOOOOD+OO 
1.900 000 OOOD+OO 
1.8000000000+00 
1.700000000D+00 
1.7900000000+00 
1 .780000000 D+00 
1.770000000D+00 
1.7600000000+00 
1.7500000 00 D+00 
1.740000000 D+00 
1.7300000000+00 
1.73 90 00 000 D+00 
1.73 8000000 D+00 
1. 73700000 OD +00 
1.7 360 00 000 D+00 
1.735000000D+00 
1 .734000000 D+00 
1.734900000 D+00 
1.734800000D+00 
1.734700000D+00 
1.7346000000+00 
1.7 34 500000 D+00 
1.734400 000 D+00 
1.734300000D+00 
1. 7342000000+00 
1.734290000D+00 


FUNCTION  D ( B * ) 

2.887203645D+00 
2 .872928 89 2 D+00 
2.854500246D+00 
2.829967336D+00 
2 .79608482 8 D+00 
2. 747210 76  ID +00 
2 .6729 59 541 D+00 
2 .552119579 D+00 
2.32972582  3D +00 
1.753356611D+00 
2 .29 6 69 5 350 D+00 
2.260528724D+00 
2.2207110950+00 
2 .2567207430+00 
2.252875780D+00 
2 .24 8993245 D+00 
2. 2450725370+00 
2.241113038D+00 
2.23711411 7D+00 
2. 23307 513 1D+ 00 
2.236712032D+00 
2.236309546D+00 
2.235906658D+00 
2. 23550336 7D+00 
2 .235099673 D+00 
2 .234 69 5 575 D+00 
2.235059282D+00 
2.235018886D+00 
2.2349784860+00 
2 .23493808  3D +00 
2 .2 34 8 9 76 7 5 D+00 
2.234857263D+00 
2.2348168470+00 
2.2347764270+00 
2.2348128050+00 


A  -  D( 8  •  ) 

-6 .5239886000-01 
-6.3812410700-01 
-6 . 196954611D-01 
-5.951625507D-01 
-5.612800428D-01 
-5.124059755D-01 
-4. 381547563 D-01 
-3. 173147936D-01 
-9 .4921037620-0 2 
4 .8144  817410-01 
-6 .1890564610-02 
-2.572393895D-02 
1 . 40936  8974 D-0 2 
-2.1915957890-02 
-1 .8070994820-02 
-1.41 88 4 60 36 D-0 2 
-1 .0267751780-02 
-6 .30 8 2 526 3 8 D-0 3 
-2 . 309332347D-0  3 
1.729654209D-03 
-1 .90 724 74 4 6 D-0  3 
-1.50476 12 41D-0 3 
-1.1018730780-03 
-6 .9858230110-04 
-2.948882514 D-04 
1. 0920972990-04 
-2.544966411D-04 
-2 .141009908  D-04 
-1 • 73701 299 9D-04 
-1.3  32975676 D-0 4 
-9. 2889793460-05 
-5. 24779766 7 D-0 5 
-1 .2062 11 661 D-0 5 
2.8357787400-05 
-8, 020 30 8 19 4 D-0 6 


ESTIMATED  LAMBDA 

1.4233760 79 D- 13 
2. 81837 44 80 D- 12 
5 .5593760180-11 
1.09056  82  8 8 D-0 9 
2.121781412D-08 
4. 0764 3 31 21 D-0 7 
7.6  8 00 000 00 D-0 6 
1.403508772D-04 
2 .4489795920-03 
4 .0000000000-02 
3. 24 8 8 791 62 D-0 3 
4.307107646D-03 
5 .7060093080-03 
4.430107818D-03 
4.556588426D-03 
4. 6866 46926 D- 03 
4 .8203834760-03 
4 .9579010000-03 
5.099 30 52720-03 
5 .2 4470498 8 D-0 3 
5.113663630D-03 
5.128062049D-03 
5 .1425006420-03 
5 .1569795180-03 
5.171498790D-03 
5 .186058568  0-03 
5. 172 952943 D-03 
5 .1744075010-03 
5.1758624650-03 
5. 1773178340-03 
5. 178773609 D-03 
5 .1802297890-03 
5 .1816863750-03 
5.183143366D-03 
5.1818320 56D-03 
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THE  FINAL  ESTIMATE  OF  BETA  -  1 ,7342900D«-00 

THE  FINAL  ESTIMATE  OF  LAMBDA  -  5 . 18 18 32 ID-03 

CONVERGENCE  TO  8 . 0203082D-06 

WHICH  IS  LESS  THAN  EPSILON  *  1 .OOOOOOOD-O 5 

THE  FINAL  STEP  SIZE  IS  1 .OOOOOOOD-05 


48 


CRAMER  -  VON  HISES  GOODNESS  OF  FIT  TEST. 

ESTIMATED  BETA  «  1 .7342900D+00 

NUMBER  OF  FAILURES  -  3 

UNBIASED  ESTIMATE  OF  BETA  «  8 .76123780-01 


CRAMER  -  VON  MISES  STATISTIC  »  5. 43057580-02 

REJECT  THE  WEIBUll  INTENSITY  MODEL  IF  THE  STATISTIC 
EXCEEDS  THE  APPROPRIATE  CRITICAL  VALUE  FOR  M  »  3 
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DATA  INPUT  PHASE 


TOTAL  NUMBER  OF  FAILURES  -  5 

TOTAL  NUMBER  OF  SYSTEMS  »  I 

SYSTEM  STARTING  AND  ENDING  TIMES. 


FAILURES. 


l.OOOD+OO  1.000D+02 


3.200D+01  A.AOOD+Ol  5.600D+01  7.500D+01  9.500D+01 


ESTIMATION  OF  THE  PARAMETERS  OF  THE  WEI8ULL  INTENSITY  FUNCTION 


CONSTANT  NOT  INVOLVING  8ET A*  A  •  A *02932840+00 


ESTIMATED  8ET A 

1.000000000D+01 
9 • 000 00 0000 D +00 
8*0000000000+00 
7.0000000000+00 
6*0000000000+00 
5  *0000000000+00 
4*0000000000+00 
3*0000000000+00 
2.0000000000+00 
l.OOOOOCOOOD+OO 
1.9000000000+00 
1*8000000000+00 
1*7000000000+00 
1.7900000000+00 
1. 7800000000+00 
1*  7 700 00000 D  +  00 
1* 7600000 00 D+GO 
1.7500000000+00 
1*7400000000+00 
1.7300000000+00 
1.7390000000+00 
1.7380000000+00 
1.7370000000+00 
1.7360000000+00 
1.7350000000+00 
1.7340000000+00 
1.7330000000+00 
1.7320000000+00 
1.7310000000+00 
1.7319000000+00 
1.731 800000 D +00 


FUNCTION  DOM 

4 .5051701860+00 
4  .4940590750+00 
4.4801701860+00 
4.4623130430+00 
4.4385035190+00 
4.4051701860+00 
4.3551702320+00 
4.2718414580+00 
4.1056307490+00 
3.6516870570+00 
4.079584382D+00 
4.0507716880+00 
4.0187689730+00 
4.0477225700+00 
4.0446411520+00 
4.0415269630+00 
4.0383795270+00 
4.0351963580+00 
4.0319829630+00 
4 .028  732840  D +00 
4.0316595210+00 
4.0313357320+00 
4.0310115950+00 
4.0306871090+00 
4 .0303622740+00 
4.0300370890+00 
4.0297115530+00 
4.0293856670+00 
4.0290594290+00 
4.0293530590+00 
4.0293204480+00 


A  -  D  C  B  •  ) 

-4.758417395D-01 
-4 . 647306284  D-0 1 
-4.50  8417395  0-0 1 
-4.3298459660-01 
-4.0917507280-01 
-3. 7584 17399 D-01 
-3.258417855 D- 01 
-2.4251301130-01 
-7.6302302550-02 
3.7764139000-01 
-5 .025 593 598 D-0 2 
-2.1443241000-02 
1.0559473210-02 
-1.8394123450-02 
-1.5312705290-02 
-1.219851680D-02 
-9.0510805740-03 
-5.8699113840-03 
-2.654516074D-03 
5.9560654950-04 
-2.3310748810-03 
-2.0072859160-03 
-1 .683148672 C- 03 
— 1.358662642 D- 03 
-1.0338273180-03 
-7.0864219200-04 
-3.8310675290-04 
— 5. 7220 49 074 D- 05 
2.690171060  D-0 4 
-2.4612549570-05 
7 . 998905460  D-0 6 


ESTIMATEO  LAM0OA 

5 .0000000000-20 
5. 000000000 D- 18 
5.0000000000-16 
5 .000000000 D- 14 
5.0000000000-12 
5 .000000001 D-10 
5 .0000000500-08 
5. 000005000 C-06 
5. 0005000500-04 
5 .05050 50 51D -02 
7. 925722105 D- 04 
1.2562587740-03 
1. 9913286150-03 
1.31548 00 02D-03 
1.377493745D-03 
1.442431753D-03 
1.5104320030-03 
1.5816389880-03 
1.65620 40 28D-03 
1 .7342 85 5 93 D-0 3 
1.6638512620-03 
1.6715338180-03 
1.6792518590-03 
1.687005549D-03 
1.6947950520-03 
1.7026205340-03 
1.7104821620-03 
1. 718380103 D- 03 
1.726314523D-03 
1.7191719000-03 
1.7199640630-03 
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THE  FINAL  ESTIMATE  OF  BETA  -  1. 73180000+00 

THE  FINAL  ESTIMATE  OF  LAMBDA  -  1 . 71996 AID-03 

CONVERGENCE  TO  7.9989055D-06 

WHICH  IS  LESS  THAN  EPSILON  ■  1.0000000D-05 

THE  FINAL  STEP  SIZE  IS  1. 0000000D-04 


52 


CRAMER  -  VON  MISES  GOODNESS  OF  FIT  TEST 


ESTIMATED  BETA  -  1.7318000D+00 

NUMBER  OF  FAILURES  •  5 

UNBIASED  ESTIMATE  OF  BETA  -  1. 3892706D+00 


THE  CRAMER  -  VON  MISES  GOODNESS  OF  FIT  TEST  IS  NOT  APPROPRIATE  FOR  THIS  CASE 
BECAUSE  ONE  OR  MORE  SYSTEMS  HAVE  NON-ZERO  STARTING  TIMES. 


DATA  INPUT  PHASE 


INDATA  ERROR 
DOES  NOT  FALL 

--  THE  FAILURE  AT  1 .OOOD+Ol 
WITHIN  THE  TEST  PERIOD  FROM 

1 .2200+01 

TO 

4. 530D+01 

INDATA  ERROR 
DOES  NOT  FALL 

—  THE  FAILURE  AT  4.579D+02 
WITHIN  THE  TEST  PERIOD  FROM 

1. 2200+01 

TO 

4 .5300+01 

TOTAL  NUMBER  OF  FAILURES  ■  A 

TOTAL  NUMBER  OF  SYSTEMS  ■  1 


SYSTEM  STARTING  AND  ENDING  TIMES. 


FAILURES. 


1.220D+01  4.530D+01 


1  .OOOD+Ol  1.450D+01  1.60OD+O1  4.579D+02 
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DATA  INPUT  PHASE. 


TOTAL  NUMBER  OF  FAILURES  «  AO 
TOTAL  NUMBER  OF  SYSTEMS  »  1 


SYSTEM  STARTING  AND  ENDING  TIMES. 


FAILURES. 


0. 


3.256D+03 


7.000D-01 
9.920D+01 
1. 745D+02 
4.905D+02 
7 • 859D+02 
1. 1360+03 
1.572D+03 
2 . 525D+03 


3.700D+00 
1. 122D+02 
1. 916D+02 
5.133D+02 
8 • 870D+02 
1. 1790+03 
1.630D+03 
2 • 929D+03 


1.3200+01 
1.209D+02 
2. 8280+02 
5. 5840+02 
1.0110+03 
1.2600+03 
1.7020+03 
3.0160+03 


1.7600+01 

1.5100+02 

3.5520+02 

6.781D+02 

1.0290+03 

1.2980+03 

1.9290+03 

3.1810+03 


5.4500+01 

1.6300+02 

4.8630+02 

6.8800+02 

1.0340+03 

1.4200+03 

2.0720+03 

3.2560+03 
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ESTIMATION  OF  THE  PARAMETERS  OF  THE  WE  I  BULL  INTENSITY  FUNCTION 


CONSTANT  NOT  INVOLVING  SETA  J  A  -  6 .04649880+00 


ESTIMATED  SETA 

l.OOOOOOOOOD+Ol 
9  <000  00  000004-00 
8. 000 000000 D+00 
7.00000000004-00 
6.00000000004-00 
5.00000000004-00 
4.000000000D+00 
3.00000000004-00 
2 .000000 000 D 4-00 
1.000000000D4-00 
9. 000000000 D-01 
8. 000 000000 D-01 
7. 000000000 D-01 
6 .000000 000 D-01 
5.0000000000-01 
A. 0 00000000 D-01 
A. 9 00 000 00 00-01 
A.  8000000000-01 
4 .8 90000 000 D-01 
A.  8990000000-01 
A ,898000000 D-01 
A. 897000000 D-01 
A. 8979000000-01 
A .8 97 800 000 D-01 
A.  897 70 00000-01 
A.  897600000 D-01 
A • 8 97 5 00 000 D-01 


FUNCTION  0 ( B  •  ) 

7.9883468600+00 
7.977235749D+00 
7.963346860D+00 
7.945489718D+00 
7. 9216801940+00 
7.888346860D+00 
7. 8383468  60  D4-00 
7.75501352704-00 
7.588346860D4-00 
7. 088346  86  OD  4-00 
6.977235749D4-00 
6.83834686004-00 
6.65977543204-00 
6.421680194D4-00 
6  .08834686004-00 
5.588346860D4-00 
6.047530534D4-00 
6  .005013527D4-00 
6.043357085D4-00 
6.047113956D4-00 
6.046697208D4-00 
6.046280289D4-00 
6.04665552304-00 
6  .04661 383 7D4-00 
6.046572150D4-00 
6. 046  530461 D4-00 
6.046488770D4-00 


A  -  D(  B  •  ) 

-1.94184810104-00 
-1.93073699004-00 
-1.91684810104-00 
-1.89899095804-00 
-1.87518143404-00 
-1.84184810104-00 
-1.79184810104-00 
-1.70851476804-00 
-1.541848101D4-00 
-1.04184810104-00 
-9,3073698980-01 
-7.9184810090-01 
-6.1327667230-01 
-3.7518143420-01 
-4.1848100910-02 
4.5815189910-01 
-1. 031774382 D- 03 
4.1485232420-02 
3.1416741390-03 
-6.151962382D-04 
-1. 984479930 D-04 
2. 184704576D-04 
-1.5676380890-04 
-1.1507792270-04 
-7.3390334170-05 
-3.1701043280-05 
9.9899500850-06 


ESTIMATEO  LAMBDA 

2.9841173600-34 
9.7171813580-31 
3.1642057660-27 
1.0303603230-23 
3.3551623210-20 
1.0925415070-16 
3.5576429080-13 
1.1584752600-09 
3.7723429900-06 
1.2283880480-02 
2.7580873920-02 
6.1927060240-02 
1.3904420870-01 
3.1219457050-01 
7.0096734530-01 
1. 5738749660+00 
7.6002001730-01 
8.2404755450-01 
7.6619225070-01 
7.6063499650-01 
7.6125047340-01 
7.6186644820-01 
7.6131204840-01 
7.6137362850-01 
7.6143521350-01 
7.6149680350-01 
7.6155839850-01 
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THE  FINAL  ESTIMATE  OF  BETA 


4 . 8975000D— 01 


THE  FINAL  ESTIMATE  OF  LAMBDA  -  7.61558400-01 

CONVERGENCE  TO  9.9899501D-06 

WHICH  IS  LESS  THAN  EPSILON  -  1 .00000000-05 

THE  FINAL  STEP  SIZE  IS  1. OOOOOOOD-05 
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CRAMER  -  VON  MISES  GOODNESS  OF  FIT  TEST 


ESTIMATED  BETA  ■  4.8975000D-01 

NUMBER  OF  FAILURES  -  40 

UNBIASED  ESTIMATE  OF  BETA  ■  4 .6526478D-01 


CRAMER  -  VON  MISES  STATISTIC  -  6. 8282282D-02 

REJECT  THE  WEIBULL  INTENSITY  MODEL  IF  THE  STATISTIC 
EXCEEDS  THE  APPROPRIATE  CRITICAL  VALUE  FOR  M  -  39 
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DATA  INPUT  PHASE 


INDATA  ERROR  —  A  FAILURE  AT  0.0000000  WAS  INPUT. 

THE  PROBABILITY  OF  SUCH  A  FAILURE  TIME  IS  0.0  ACCORDING  TO  THE  MODEL. 

TOTAL  NUMBER  OF  FAILURES  -  7 

TOTAL  NUMBER  OF  SYSTEMS  »  1 

SYSTEM  STARTING  AND  ENDING  TIMES. 

FAILURES. 

0.  7. 9000+01 

0.  2.0000+01  4.500D+01  6.2000+01  6.4000+01 

7 .9000+01  6.8000+01 
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DATA  INPUT  PHASE 


TOTAL  NUMBER  OF  FAILURES  - 

5 

TOTAL  NUMBER  OF  SYSTEMS  - 

1 

SYSTEM  STARTING  AND  ENDING 

TIMES • 

FAILURES. 

4.560D+02  7 • 890  D+02 

4.560D402  4.561D+02  4.570D402  4.580D402  4.fi90D+02 

* 
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ESTIMATION  OF  THE  PARAMETERS  OF  THE  WEI  BULL  INTENSITY  FUNCTION 


CONSTANT  NOT  INVOLVING  BETA*  A  -  6.13782400+00 


ESTIMATED  BETA 

1.0000000000+01 
9.000000000D+00 
8.0000000000+00 
7.0000000000+00 
6.0000000000+00 
5.0000000000+00 
4.0000000000+00 
3.0000000000+00 
2.000000000D+00 
1.0000000000+00 
9.0000000000-01 
8.0000000000-01 
7.0000000000-01 
6.0000000000-01 
5.0000000000-01 
4 .000000000 D-01 
3.0000000000-01 
2.0000000000-01 
1.0000000000-01 
9.0000000000-02 
8.0000000000-02 
7.0000000000-02 
6.0000000000-02 
5.0000000000-02 
4.0000000000-02 
3.0000000000-02 
2.0000000000-02 
1.0000000000-02 
9.0000000000-03 
8.0000000000-03 
7.0000000000-03 
6.0000000000-03 
5.0000000000-03 
4. 0000000000-03 
3.0000000000-03 
2.000000000D-03 
1.0000000000-03 
9.0000000000-04 
8.0000000000-04 
7.0000000000-04 
6.0000000000-04 
5.0000000000-04 
4.0000000000-04 
3.0000000000-04 
2.0000000000-04 
1.0000000000-04 
9.0000000000-05 
8.0000000000-05 
7.0000000000-05 
6.0000000000-05 


FUNCTION  0 ( B • ) 

6.5730555290+00 
6.5636282560+00 
6.5526773070+00 
6.5399780970+00 
6.5253232290+00 
6.5085569710+00 
6.4896198020+00 
6.4685963840+00 
6.4457540820+00 
6.4215552730+00 
6.4190838880+00 
6.4166058560+00 
6.4141218920+00 
6.4116327180+00 
6.4091390650+00 
6.4066416700+00 
6.4041412750+00 
6.4016386250+00 
6.3991344720+00 
6.3988  8  400  30  +00 
6.398633527D+00 
6.3983830450+00 
6.3981325570+00 
6.3978820660+00 
6.3976315700+00 
6.3973810710+00 
6.3971305710+00 
6.3968800680+00 
6.3968550180+00 
6.3968299680+00 
6.3968049170+00 
6.3967798670+00 
6.3967548170+00 
6.396729766D+00 
6.3967047160+00 
6.3966796660+00 
6.3966546160+00 
6.3966521100+00 
6.3966496050+00 
6.3966471000+00 
6.3966445950+00 
6.3966420900+00 
6.396639585D+00 
6.3966370800+00 
6. 396634575 D  +  00 
6.3966320700+00 
6.3966318200+00 
6.3966315690+00 
6.3966313190+00 
6.3966310680+00 


A  -  0  (  B  •  ) 

-4.3523153700-01 
-4.2580426400-01 
-4 .148533158 D— 01 
-4.0215410580-01 
-3.8749923  7  8  D-01 
—3. 70732  9  79  30—0 1 
-3.5179581000-01 
-3. 3077239200-01 
-3.0793009020-01 
-2.8373128180-01 
-2.8125989690-01 
-2.7878186490-01 
-2.7629790050-01 
-2.7380872670-01 
-2.7131507380-01 
-2.6881767870-01 
-2.6631728330-01 
-2.6381463390-01 
-2.6131048020-01 
-2.6106001100-01 
-2.6080953500-01 
-2. 6055905300-01 
-2.6030856570-01 
-2.6005807400-01 
-2.5980757840-01 
-2.5955707990-01 
-2.5930657900-01 
-2.5905607670-01 
-2.5903102640-01 
-2.5900597610-01 
-2.589 80 92 58D -01 
-2.5895587550-01 
-2.5893082520-01 
-2.5890577490-01 
-2.5888072460-01 
-2.5885567430-01 
-2.5883062400-01 
-2.5882811890-01 
-2.5882561390-01 
-2.5882310890-01 
-2.5882060380-01 
-2.5881809880-01 
-2.5881559380-01 
-2.5881308870-01 
-2.5881058370-01 
-2. 5880807870-01 
-2.5880782820-01 
-2.5880757770-01 
-2.5880732720-01 
-2.5880707670-01 


ESTIMATEO  LAMBOA 

5.3704355370-29 
4.2502328940-26 
3.3712739540-23 
2.6846472360-20 
2.1527927750-17 
1.7479656280-14 
1  .4522429250-11 
1.2615104850-08 
1.2060253020-05 
1.5015015020-02 
3.1704011970-02 
6.7762626330-02 
1 .4  70949998  D-0 1 
3.2587604730-01 
7.4239189660-01 
1.7613042810+00 
4.4561084710+00 
1. 2680050910+01 
4. 8096807 79 D+01 
5.6972370190+01 
6.8329189010+01 
8.325 04 4844D+01 
1.0354297810+02 
1.3246101870+02 
1.765156468D+02 
2.5090306230+02 
4.0121776100+02 
8.5544479230+02 
9.5659388780+02 
1.0830742670+03 
1.2457425170+03 
1.4626929500+03 
1.7664953360+03 
2.2222892110+03 
2.9820668400+03 
4.5018050160+03 
9.0613877230+03 
1.0074650930+04 
1.1341234560+04 
1. 29697045 4D+04 
1.5141004030+04 
1.8180830750+04 
2.2740580110+04 
3.0340174790+04 
4.5539382760+04 
9.1137043900+04 
1.0126985980+05 
1.1393588010+05 
1.3022076380+05 
1.5193394280+05 
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5. 000000000 D- 05 
4. 000 000000 D-0 5 
3.000000000D-05 
2. 000 000000 D-05 
1 *000 000000 D-0 5 
9.000000000D-06 
8.000000000D-06 
7. 000000 000 D-0 6 
6 >0000 00000 D-0 6 
5. 000 0000 00 D- 06 
4. 0000000 00D- 06 
3  *000000000 D-0 6 
2 *000000000 D-0 6 
1* 00 0000000 D-06 
9, 000 000000 D-0 7 
8. 000 000 000 D-0 7 
7 • 000000 000 D-0 7 
6.000000000D-07 
5.000000000D-07 
4.000000000D-07 
3.000000000D-07 
2. 000000000 D-07 
1  .OOOOOOOOOD-07 
9 .00000 0000 D-0 8 
8.000000000D-08 
7. 000 000000 D- 08 
6. 000000 000 D-0 8 
5.000000000D-08 
4. 000 000 000 D— 08 
3.000000000D-08 
2.000000000D-08 
1.000000000D-08 


6.396630818D+00 
6.396630567D+00 
6.396630317D+00 
6 .3 9663 006 6D+ 00 
6.396629816D+00 
6. 396629 79  ID ♦00 
6.396629766D+00 
6.396629741D+00 
6.396629715D+00 
6.396629690D+00 
6.396629665D+00 
6. 396629 64 OD +00 
6.396629615D+00 
6.396629590D+00 
6, 396629 58 8D +00 
6.396629585D+00 
6.396629583D+00 
6.396629580D+00 
6.396629578D+00 
6.396629575D+00 
6. 39662957 3D +00 
6.3  96629570D+00 
6. 39662956 8D+00 
6 .39662956  7D+00 
6.396629567D+00 
6.3  9662 9 56 7D +00 
6 . 3 9662 9 56 7D +00 
6. 396 62 9 56 6D +00 
6 . 396629  566 D+00 
6. 396 62 956 6D +00 
6.396629566D+00 
6.396629  565D+00 


-2. 588068262 D-01 
-2, 588065756D-01 
-2. 5 8 806 3 251 D-01 
-2.588060746D-01 
-2. 588058241 D-01 
-2 . 5  8805  799 1 D-01 
-2.588057740D-01 
-2. 5 88057490 D-01 
-2.588057239D-01 
-2. 588056989 D-01 
-2. 5 88056738 D-01 
-2.588056488D-01 
-2.588056237D-01 
-2.58 8055987 D-01 
-2.588055962D-01 
-2.588055  9  3  7D-01 
-2 . 588055912 D-01 
-2, 588055887D-01 
-2 .58 80 55 862 D-01 
-2. 588055837 D-01 
-2. 5 8 80 55 812 D-01 
-2.588055786D-01 
-2.588055761D-01 
-2, 5 8 80 5 5 759 D-01 
-2. 5 8 805 5 75 6 D-01 
-2.588055754D-01 
—2. 5 8 80 5 5751 D-01 
-2, 5 88 05 5 749 D-01 
-2. 588055746D-01 
-2.588055744D-01 
-2.588055741D-01 
-2.5  880 5 5 73 9 D-01 


ML EW  ERROR  --  BETA  LESS  THAN  0.00000001 


1.823323941D+05 
2.279300721D+05 
3. 03926 20 31D +05 
4.559184672D+05 
9.118952631D+05 
1.0132234 40 D+06 
1.139883662D+06 
1.302732518D+06 
1.519864326D+06 
1.823848858D+06 
2.279825655D+06 
3.039786985D+06 
4.559709644D+06 
9.119477621D+06 
1.013275939D+07 
1.139936161D+07 
1.302785017D+07 
1.519916826D+07 
1.823901358D+07 
2.279878155D+07 
3.039839485D+07 
4.559762144D+07 
9.119530122D+07 
1.0132811 89 D+0 8 
1.139941411D+08 
1.302790267D+08 
1.519922076D+08 
1.823906608D+08 
2.279883406D+08 
3.039844735D+08 
4.559767394D+08 
9.119535372D+08 


62 


DATA  INPUT  PHASE 


TOTAL  NUMBER  OF  FAILURES  »  12 

TOTAL  NUMBER  OF  SYSTEMS  »  1 


SYSTEM  STARTING  AND  ENDING  TIMES. 


FAILURES. 


8  .QOOD  +  02 

4  •500D+01 

4*56004*02 

A ,6700+02 

4.77004*02 

4.64004*02 

4.920D4-02 

1.2000-01 

5. 6400+02 

6.50004*01 

7.60004*01 

8  *  900D+01 

9.91004*01 
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ESTIMATION  OF  THE  PARAMETERS  OF  THE  WEIBUll  INTENSITY  FUNCTION. 
CONSTANT  NOT  INVOLVING  SET  A  J  A  -  4.70451810+00 


ESTIMATED  BETA 

l.OOOOOOOOOD+Ol 
9. OOOOOGOOOD+OO 
8. OOO OOOO 00 D +00 
7 .0000000000+00 
6.000000000D+00 
5.000000000D+00 
4.000000000D+00 
3.000000000D+00 
2. 0000000000+00 
l.OOOOOOOOOD+OO 
9 • 000000 000 D— 01 
8.000000000D-01 
7,  0000000000-01 
6.0000000000-01 
5.000000000D-01 
5.900000000D-01 
5.8000000000—01 
5.7000000000-01 
5. 6000000 00 D-01 
5.5000000000-01 
5. 400000000 D-01 
5.3000000000-01 
5 . 2000 000 00 D— 01 
5.1000000000-01 
5.000000000D-01 
5 ,090 000000 D-01 
5.0800000000-01 
5  .070000 000D-01 
5.0600000000-01 
5  .050000000D-01 
5 .059 000000 D-01 
5.0580000000-01 
5.0570000000-01 
5.0560000000-01 
5.055000000D-01 
5, 054000000 D-01 
5.0530000000-01 
5  .0520 00 000 D-01 
5.0510000000-01 
5. 0500000 00 D-01 
5.0509000000-01 
5.0508000000-01 
5.0507000000-01 
5.050600000D-01 
5  .050  500000  D-0 1 
5 • 050400000D-01 
5.0503000000-01 
5.0502000000-01 
5.0502900000-01 


FUNCTION  0 ( B • ) 

6.584611728D+00 

6 • 57350061 7 D+00 

6 , 5596 1172 8 D +00 

6.5417545850+00 

6.5179450610+00 

6.484611728D+00 

6, 434611 72 8 D+00 

6.3512783940+00 

6.1846117280+00 

5,6846117280+00 

5, 573 500617 D+00 

5.434611728D+00 

5.256040299D+00 

5.0179450610+00 

4, 6 846 1172 8 D+00 

4. 989696473 D+00 

4.960473797D+00 

4 .930225763D+00 

4.8988974420+00 

4.8664299090+00 

4. 8 32 7 59 876 D+00 

4.7978192750+00 

4. 76 15 3480 5D +00 

4.723827414D+00 

4.684611728D+00 

4. 71997 518 5 D+00 

4.716107791D+00 

4.712225140D+00 

4.708327143D+00 

4.7044137080+00 

4 .7079364950+00 

4.7075456940+00 

4.707154737D+00 

4.7067636260+00 

4 . 706372361 D+00 

4.7059809400+00 

4.7055893650+00 

4 .705197634D+00 

4.7048057490+00 

4 .7044137080+00 

4.704766552D+00 

4 .7047273530+00 

4.7046881530+00 

4. 70464895 1D+00 

4.704609748D+00 

4.704570543D+00 

4  .7045313360+00 

4.7044921280+00 

4.7045274160+00 


A  -  0  (  B '  ) 

-1.8800936200+00 
-1.8689825090+00 
-1.8550936200+00 
-1  .837236477D+00 
-1.8134269530+00 
-1.7800936200+00 
-1. 7300936200+00 
-1.6467602870+00 
-1 • 4 80 09 3 6 20 D+00 
-9.8009361980-01 
-8.6898250870-01 
-7.3009361980-01 
-5.5152219130-01 
-3. 134269532 D-01 
1.990638017D-02 
-2.8517836560-01 
-2.5595568880-01 
-2.2570765490-01 
-1.9437933410-01 
-1.619118017D-01 
-1.2824176800-01 
-9.3301167000-02 
-5.7016696760-02 
-1.9309306110-02 
1, 9 90638017 D-02 
-1.5457077590-02 
-1.158968283  D-02 
-7.707032063D-03 
-3.809034854D-03 
1.043999676D-04 
-3.4183875750-03 
-3.0275858290-03 
-2.6366295240-03 
-2.2455185690-03 
-1 ,8 54 2 52871 D-0 3 
-1.4628323390-03 
-1.0712 5688 2D-03 
-6, 7 9 52 64 061 D-0 4 
-2.8764082040-04 
1.0439996760-04 
-2.4844372720-04 
-2.0924508190-04 
-1.7004488440-04 
-1.3084313460-04 
-9.1639832430-05 
-5.2434977770-05 
-1.3228570530-05 
2.5979389380-05 
-9.3078444090-06 


ESTIMATED  LAMBDA 

1.1175870900-28 
8.9406967160-26 
7.1525573730-23 
5, 722045898 D- 20 
4.5776367190-17 
3.6621093750-14 
2.9296875000-11 
2.3437500000-08 
1.8750000000-05 
1.5000000000-02 
2.9268485990-02 
5.7109618160-02 
1.1143413730-01 
2. 174338991 D-01 
4.2426406870-01 
2.3246531000-01 
2.4853585650-01 
2.6571737500-01 
2.8408666820-01 
3.0372584790-01 
3.2472270280-01 
3.4717109000-01 
3.7117135540-01 
3  .9683078190-01 
4.2426406870-01 
3.9949232730-01 
4.021717238D-01 
4.0486909100-01 
4.0758454950-01 
4. 1031822050-01 
4.0785709500-01 
4.0812982280-01 
4.0840273290-01 
4.0867582560-01 
4.0894910080-01 
4. 092225588 D-01 
4.0949619960-01 
4.097700235D-01 
4.1004403040-01 
4.103182205D-01 
4.100714412D-01 
4.100988538D— 01 
4.1012626820-01 
4.101536845D-01 
4.1018110260-01 
4.102085225D-01 
4. 10235 94 42D -01 
4.1026336780-01 
4.102386865D-01 
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THE  FINAL  ESTIMATE  OF  BETA  -  5 . 0502900D-01 

THE  FINAL  ESTIMATE  OF  LAMBDA  -  A .1023869D-01 

CONVERGENCE  TO  9. 3078444D-06 

WHICH  IS  LESS  THAN  EPSILON  -  1 .OOOOOOOD-05 

THE  FINAL  STEP  SIZE  IS  1 . OOOOOOOD-06 
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CRAMER  -  VON  MISES  GOODNESS  OF  FIT  TEST. 

ESTIMATED  BETA  ■  5 .0502900D-01 

NUMBER  OF  FAILURES  -  12 

UNBIASED  ESTIMATE  OF  BETA  -  4.6294107D-01 


CRAMER  -  VON  MISES  STATISTIC  -  1 .4645334D-01 

REJECT  THE  WEIBUll  INTENSITY  MODEL  IF  THE  STATISTIC 
EXCEEDS  THE  APPROPRIATE  CRITICAL  VALUE  FOR  M  ■  12 
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DATA  INPUT  PHASE 


TOTAL  NUMBER  OF  FAILURES  - 

13 

TOTAL  NUMBER  OF  SYSTEMS  - 

A 

SYSTEM  STARTING  AND  ENDING 

TIMES. 

FAILURES. 

0.  1*0000401 

0.  3.0000401 

2.0000400  3.0000400  4.0000400  5.0000400  6.0000400 
8.0000400 

0.  5.0000401 

NO  FAILURES 

0.  6.000D401 

4.1000401  4.2000401  4.3000401  4.4000401 

5.1000401  5.2000401  5.3000401 
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ESTIMATION  OF  THE  PARAMETERS  OF  THE  WEI  BULL  INTENSITY  FUNCTION 


CONSTANT  NOT  INVOLVING  BET A«  A  -  2 . 73 143730*00 


ESTIMATED  BETA 

1.000000000D+01 
9.000000000D+00 
8. 000 000000 D+ 00 
7.000000000D+00 
6.000000000D+00 
5*0000000000+00 
4.0000000000+00 
3. 000000 000 D+00 
2*0000000000+00 
1.0000000000+00 
9, 0000 00000 D -01 
9. 900000000 D-01 
9.800000000D-01 
9. 700000000D-0 1 
9.600000000D-01 
9.690000000D-01 
9.680000000D-01 
9.670000000D-01 
9.660000000  D-0  1 
9.650000000D-01 
9. 6 590000 00 D-01 
9.658000000D-01 
9.657000000D-01 
9. 65600000 0D-01 
9.655000000D-01 
9. 6 54 000000 D-01 
9.653000000D-01 
9. 652 000 000 D—0 1 
9.651000000D-01 
9. 6 50 00 00 00 D-0 1 
9.650900000D-01 
9, 6 50 80 0 000 D— 01 
9. 6 50 7000 00 D-01 
9. 6 50 60 0000 D-01 
9. 650500000 D-01 
9.650400000D-01 
9. 6 50 3000 00 D-01 


FUNCTION  D ( B * ) 

3*96843 20 09D+00 
3  .952550770D+00 
3.932861136D+00 
3.90773552 5D+00 
3.874419411D+00 
3. 8 2 794 8 74 4 D +00 
3.758529458D+00 
3.643675381D+00 
3.417047055D+00 
2 .7754906430+00 
2.641348518D+00 
2.763183818D+00 
2.7 50 649 70 4 D+00 
2.737881746D+00 
2 .7248731270+00 
2 .736591 829 D+00 
2 .735299498D  +00 
2. 7340 04 74 7 D+00 
2.732707569D+00 
2. 731407 95 6D +00 
2.732577717D+00 
2.732447B41D+00 
2.732317941D+00 
2.732188016D+00 
2 .73205 80 67D +00 
2 .7319280940+00 
2 .7317980960+00 
2.731668074D+00 
2.73153 8027D +00 
2. 73140 79 56D +00 
2 .73152502 1D+00 
2 . 731 51 201 5 D+00 
2 .7314990080+00 
2 .731486002D+00 
2 .731472995D+00 
2.731459987D+00 
2 .7 31 4469 80 D+00 


A  -  DIB1) 

-1.2369946670+00 
-1. 2211134280+00 
-1. 2014237940+00 
-1.176298182D+00 
-1. 1429820680+00 
-1.096511402D+00 
-1. 0270921150+00 
-9. 1223803830-01 
-6. 8 560 9712 4 D-01 
-4. 4053300 01 D-02 
9.0088825000-02 
-3. 1746475 79 D-02 
-1.921236129 D-02 
-6. 4 4440 36 20 D-0 3 
6.564215663D— 03 
-5.1 544  86343D-03 
-3.862155711D-03 
-2. 56 740 4 72 8 D-0 3 
-1.2702263 70 D- 03 
2.938641763D-05 
-1.140374 7570-03 
-1. 0104987930-03 
— 8.805984704D-04 
-7.5067378240-04 
-6. 2072472180-04 
— 4.907512815D-04 
-3.607534  545  D-0 4 
-2. 3073123370-04 
-1 .0068461210-04 
2. 9386417630-05 
-8. 767860 75 6D -05 
-7. 4672358980-05 
-6. 1665866310-05 
-4 .8659129540-0 5 
-3.5652148660-05 
-2.2644923660-05 
-9.6374545390-06 


ESTIMATED  LAMBDA 

1.8494583310-17 
1. 0787929630-15 
6.2596211580-14 
3. 6086151520-12 
2. 0631318340-10 
1. 1664423510-08 
6.49026 46 03D-07 
3. 5230352300-05 
1.830985915D-03 
8. 6666666670-02 
1. 2627790440-01 
9.0000321020-02 
9.3460133760-02 
9. 7050777 75 D-02 
1.0077709410-01 
9.7417222760-02 
9.7785029290-02 
9.8154202290-02 
9. 852 474672 D-02 
9.8896667580-02 
9.8561876780-02 
9.8599020620-02 
9.8636178220-02 
9.8673349610-02 
9.8710534780-02 
9.8747733740-02 
9.8784946500-02 
9.8822173050-02 
9.8859413410-02 
9.8896667580-02 
9.8863138200-02 
9.8866863140-02 
9.887058821D-02 
9.8874313420-02 
9.8878038770-02 
9.8881764250-02 
9.8885489880-02 
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THE  FINAL  ESTIMATE  QF  BETA  -  9.65030000-01 

THE  FINAL  ESTIMATE  OF  LAMBDA  -  9.8685A90D-02 

CONVERGENCE  TO  9. 637A5 A5D-06 

WHICH  IS  LESS  THAN  EPSILON  -  1 .0000000D-05 

THE  FINAL  STEP  SIZE  IS  1 . 0000000D-05 
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CRAMER  -  VON  MISES  GOODNESS  OF  FIT  TEST 


ESTIMATED  BETA  ■  9 .6503000D-01 

NUMBER  OF  FAILURES  -  13 

UNBIASED  ESTIMATE  OF  BETA  -  1.9237087D+00 


CRAMER  -  VON  MISES  STATISTIC  -  1.6179976D-01 

REJECT  THE  WEIBULL  INTENSITY  MODEL  IF  THE  STATISTIC 
EXCEEDS  THE  APPROPRIATE  CRITICAL  VALUE  FOR  M  -  13 
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OATA  INPUT  PHASE 


PROGRAM  RUN  ENOS  NORMALLY. 


Next  page  is  blank. 
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